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STUDY on using a digital ride quality augmentation 

SYSTEM TO TRIM AN ENGINE-OUT IN A CESSNA 402B 

ABSTRACT 

A linear model o£ the Cessna 402B was used to 
determine if the control power available to a Ride 
Quality Augmentation System was adequate to trim an 
engine-out. Two simulations were completed: one using a 
steady-state model, and the other using a state matrix 
model. The amount of rudder available was not 
sufficient in all cases to completely trim the airplane, 
but it was enough to give the pilot valuable reaction 
time. The system would be an added measure of safety 
for only a relatively small amount of development. 
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1. INTRODUCTION 


Due to the large increase in the commuter airline 
industry, with the federal deregulation of major 
carriers in 1978, more people are riding in small, 
short-haul, propeller-driven aircraft. This has caused 
an increased effort to make riding in such an aircraft 
as pleasant as possible. One project undertaken has 
been the development of a Ride Quality Augmentation 
System (RQAS) . This system uses acceleration feedback 
to lessen the ef'ects of turbulence by counteracting the 
undesirable accelerations with appropriate separate 
surface control deflections. 

It was felt by the developers of this system that 

since the control power was available, it would also be 

beneficial to use this system to trim the airplane in 

the event of an engine-out. This was the basis for the 

following investigation into using a RQAS to control a 

Cessna 402B in an engine-out flight condition. The 

investigation was conducted as part of NASA grant NAG1- 

345. Mr. Lou Williams is the grant monitor. 

» 

Chapter 2 of this report describes the proposed 
RQAS for the Cessna 402B. The models used to evaluate 
the system in an engine-out are then given: the steady- 
state model in Chapter 3 and the small perturbation 
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model in Chapter 4. The results of using these models 
for an engine-out are then discussed in Chapter 5. 

2. RIDE QUALITY AUGMENTATION SYSTEM 
Reference 1 documents the proposed control surface 
modifications for the RQAS in the Lateral-Directional 
mode (see Figure 1): 

1. Replace the outboard section of the split flap 
with a differentially deflecting plain flap that 
can deflect +15 to -45 degrees. 

2. Use the entire existing rudder (limiting the 
RQAS range of deflections to +5 degrees). 

The outboard flaps were computed to have 67% of the 
control power of the ailerons in the nonlinear model. 

Table 1 summarizes the control surface and actuator 
requirements for the RQAS of reference 1. 

3. LINEAR STEADY-STATE MODEL 
A nonlinear simulation model of the Cessna 402B 
(Ref. 2) was used to obtain a linear model about its 
most critical condition for an engine out. That is: 

1. sea level; 

2. flaps down; 

3. full throttle; 

4. maximum landing weight. 
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Table 1 Control Surface and Actuator 
Requirements 


Control Surface Deflection and Rate Limits 


Surface 


Deflection Maximum 

Range (deg) Rate (deg/sec) 


SS Elevator 
Flap( inboard) 
Differential Flap 
Rudder 1 


± 5 

fl5 to -45 
♦15 to -45 
♦32 


50 

120 

120 

50 


Actuator Requirements 


Surface 

Max Load 
(lbs) 

Speeed 
( in/sec) 

Stroke 

(in) 

SS Elevator 

65 

3.50 

0.75 

Flap( inboard) 

750 

8.50 

4.25 

Differential Flap 

380 

8.50 

4.25 

Rudder 

520 

3.50 

4.50 

1 The deflection 
uses a deflection 

is for the standard rudder, 
range of +5 degrees. 

The RQAS 
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This will be called 

Flight Condition #1 

1 FC1 ) and is 

defined as: 



One Engine Out 

V » 130 fps 

■ 0 ft 

Full Throttle 

X C g ■ 0.25 

W - 6200 lbs 

Full Flaps 




The stability and control derivatives given in Appendix 
A, Table A. 1 are for the maximum landing weight. 

Because the maximum landing weight differs from tH* 
maximum takeoff weight by less than 2 %, the 7 * values 
were used without correction. FI.* rri. vnditior. #1 would 
be the condition in an emerge; -y go-around. 


3.1 STEADY- STATE EQUATIONS 

The basic assumption made to determine needed 
control surface deflection was that the airplane motion 
could be modelled about a steady-state point as a set of 
first-order differential equations, as shown in Eqn. 2.1 
for an engine-out flight condition. 


' c *b c n DF C *« R 


r B 1 

f 

Cl 0 Cl « DP Cl « R 

< 

«DF 

t 

_ C «B C "« DF 


*R 1 

) 



- (mgsindcosY + Fy<r) ] 

_s 

-L«r 
qSb 


> ( 2 . 1 ) 


-(Nt ♦ ANd) 
qSg 

These equations have been uncoupled from the full 6 
degree of freedom equations by choosing bank angle, ♦. 
They are written in the stability axis system (see 


5 


section 3.2), Their derivation can be found in 
reference 3. 

3.2 COORDINATE REFERENCE SYSTEM 

The body-axis system is an orthogonal, right-hand 
set of axes with its origin at the airplane's center of 
mass. The X-axis lies along the centerline of the body. 
The X- and Z-axes lie in the airplane plane of symmetry, 
while the Y-axis is pointed out the right wing of the 
airplane. This can be seen in Figure 2. 

The non-dimensional derivatives listed in Appendix 
A, Table A.l are given in the stability-axis system. 

This system is also an orthogonal, right-hand set of 
axes with its origin at the center of mass of the 
airplane. The difference between the body- and the 
stability-axis systems is that the stability X-axis is 
oriented in the direction of the steady-state velocity 
of the airplane on its XZ-plane. They both share the 
same Y-axis as seen in Figure 2. 

3.3 DEVELOPMENT 

The thrust was calculated using an engine model and 
its average propeller efficiency. The average propeller 
efficiency is greater than the actual propeller 
efficiency by approximately 6%; therefore, the actual 
pitching moment and yawing moment during engine-out 
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Stability Axis Coordinate Syatam $ 

•t 

FIGURE 2 Airplane Coordinate Systems ■ 

J 



would be slightly less* The propeller efficiencies are 
given in Appendix A. 

By assuming a weight and flight condition, the 
angle of attack was obtained from the airplane lift- 
curve slope and intercept as shown in Appendix A. Some 
of the derivatives are functions of angle of attack. 

Once the angle of attack was calculated, the non- 
dimensional derivatives were obtained from reference 2. 
They are given in Appendix A. 

By varying the speed and the vertical tail size and 
solving Equation 2*1 as shown in Appendix A, the 
sideslip, rudder, and differential flap deflections were 
determined* Figure 3 shows how rudder deflection varies 
with flight condition. 

3.4 RESULTS 

It can be Seen from Figure 3 that at speeds below 
approximately 125 fps, there is not enough rudder to 
keep the airplane in straight-line flight. This is the 
minimum control speed, Vmc* This is 15 fps less than 
the minimum control speed' given in the operating 

a 

handbook; therefore, a minimum control speed of 130 fps 
is still conservative. At 130 fps, it was found that to 
fly straight with a bank angle, of -5° required 30° 
of rudder deflection, 7.5° of differential flap 
deflection, with a sideslip of -3°. The pedal force 
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FIGURE 3 Rudder required for en engine-out 
on a Cessna 402B 




required was also found to be within federal 
regulations. 

It can also be seen from Figure 3 that if one- 
fourth of the vertical tail were removed, the airplane 
would no longer meet FAR 23 requirements. FAR 23 
requires that the minimum control speed, Vmc be greater 
than 1.2Vg. This model does not account for any 
transient phenomena of the airplane in reaching its 
steady-state condition. 

4. SHALL PERTURBATION MODEL 
The non-linear simulation model was also used to 
develop a linear model to study the dynamic behavior of 
the Cessna 402B in open- and closed- loop simulations. 
This was only done for FC1. 

4.1 EQUATIONS OF MOTION 

From the assumptions made in the steady-state 
model, a new set of matrix equations can be written: 

x - A x + B u (3.1) 

where 

» 

x’ - {&, p, r, $}, and 

u* ■ C6pp, 6 r) . 

Derivations of this equation can be found in 
reference 4. These matrices are defined in Appendix A. 
The matrices used for FC1 are given in Table 2. 


10 


Table 2 small Perturbation Matrices 


Cessna 40 2B (PCI) 


”-0.105 

-0.000867 

-0.9899 

0.245 

-1.329 

-1.752 


0.483 

O 

• 

o 

1.408 

-0.0428 


-0.299 

o 

• 

o 

0.0 

1.0 


0.149 

0.0 J 






” 0.0 

0.030(T 




-0.968 

0.206 




0.0603 

-0.835 




0.0 

0.0 





D 


0.0 

-0.0468 

0.445 

0.0 


f 
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Equation 3,1 is written with the following assumptions: 

1. perturbations are small, and 

2. initial condition is a straight-line 
trimmed flight condition. 

4.2 COORDINATE REFERENCE SYSTEM 

All of the dimensional and non-dimensional 
derivatives were calculated in the stability-axis 
system of Figure 2. The instruments onboard the 
airplane will sense the body-axis motion. This can be 
simply transformed to the stability-axis by rotating 
about the Y-axis by the airplane angle of attack. 

4.3 DEVELOPMENT 

To simulate an engine-out situation using the state 
matrices, a disturbance matrix, D, was added to the 
state and control matrices in Equation 3.1 yielding 
Equation 3.2. This disturbance matrix was made up of 
the constant angular accelerations imparted on the 
airplane in the pitch and yaw directions due to the 
engine-out. In all the cases, this matrix was commanded 
to "turn on" at one second into the simulation. 

t 

x*Ax+Bu+Dw (3.2) 

where 

w » {0 or 1} 

The open-loop response of the airplane with the 
addition of the disturbance matrix was calculated using 
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the Interactive Control Augmentation Design program 
( I cad) [ 5] • This was done on the University of Kansas 
School of Engineering's Harris computer system. The 
closed-loop response was also done on ICAD. 

As allowed by FAR 23 , a bank angle of -5° was used 
to lower the rudder required as much as possible. This 
turned out to be a difficult task on ICAD as currently 
written. Ideally, bank angle, sideslip, and yaw rate 
should be driven to -5°, 0, and 0 simultaneously. To 
get these, only bank angle was commanded and all the 
variables except roll rate were weighted heavily. This 
gave large transient values for control positions and 
rates, but in an actual system this would not happen. 
Figure 4 shows an example time history. The average 
values of control deflections and control rates have 
been drawn to show an approximation of what would 
actually happen. 

It can be seen that to trim the airplane requires 
all 32° of the rudder to be available to the RQAS. As 
proposed, only 5° of rudder deflection is available to 

t 

the RQAS; therefore, time histories were also simulated 
allowing only 5 and 15 degrees of rudder deflection. 
With 15° of rudder, the RQAS would be able to trim an 
engine out at speeds above about 170 fps. Five degrees 
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FIGURE 4 Time history of an engine-out in a 
RQAS controlled Cessna 402B 


only slows the rates and could not trim the airplane 
below its "never exceed" speed. 

5. DISCUSSION OF RESULTS 

The control deflections obtained for FC1 were: 

Model B 6 r 6np 

(deg) (deg) (deg) 

Steady-State -2.9 30.4 7.5 

Small Perturbation -5.0 30.0 5.0 

These values were obtained for a steady-state velocity 
of 130 fps. The difference in the values between the 
two models is due to the fact that in the small 
perturbation model, yaw and pitch rate are not exactly 
zero. The bank angle is also not quite -5°. It can 
also be noted from Figure 4 that the RQAS has the 
airplane trimmed in a fraction of a second, much faster 
than a human pilot could react. This result is obtained 
assuming the servo is modelled as a first-order lag. 

The slowing of the yaw rate due to only 5° of 
rudder deflection would give the pilot valuable time to 
react to the increasing sideslip. In addition, the 
differential flaps have far more control power than is 
needed to control roll and bank angle. The quicker the 
airplane reaches a bank angle of -5°, the longer the 
pilot has to respond to the yaw. 


15 


The RQAS has two obvious advantages. First, if the 
dynamic pressure is high enough and the available rudder 
is adequate , the RQAS could control an engine-out 
without the pilot touching the controls. Second, the 
RQAS would give the pilot more time to react to the 
engine-out even if the RQAS was not able to trim the 
airplane entirely. 

In addition, the RQAS will introduce digital 
computers into commuter airplanes which could then be 
used for many other jobs which of themselves do not 
warrant the expense of a computer. 

If only 5° of rudder deflection is deemed adequate, 
no modifications need to be made in addition to those 
proposed for the RQAS. The RQAS block diagram is shown 
in Figure 5. This system would treat the engine— out as 
a large disturbance. The gains chosen in the normal use 
the RQAS, however, might not be suitable for the 
engine-out. In that c&se, engine-out sensors would be 
required to tell the system when to change gains. 

Giving more conti ol power to the system for an 
engine-out could be done in two ways. First, the RQAS 
could be given the authority to move the rudder more 
than 5*. This would require a proportionately stronger 
actuator, but it would not change the stroke given in 
Table 1 which is for 32°. Secondly, an engine-out 
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FIGURE 5 RQAS Block Diagram 
















sensor could trigger s change in the amount of rudder 
allowed so that the RQAS is still only allowed 5* while 

the engine-out is allowed mores 

No additional software will be required If gain 
scheduling is not needed or is already incorporated in 
the RQAS. The components needed for this system are 
hydraulic pump to be mounted on one of the engines, an 
accumulator, and a set of electro-hydraulic actuators. 
The design of the actuator and its implementation will 
be contracted to Cessna. In addition, a set of engine- 
out sensors will be required. Reliability, ease of 
implementation, and cost of these should be 
investigated. 

6. CONCLl f fiTftWS AMD RECOMMENDED RESEARCH 
6.1 CONCLUSIONS 

The following conclusions were reached: 

1. The RQAS with control of all 32° of rudder 
deflection can trim an engine-out in any steady- 
state flight condition. 

2. With only th. proposed 5" of rudder deflection, 
the RQAS can slow the divergence of the airplane 
significantly. 

3. Using only the proposed control surface 
modifications, implementing the engine-out 
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capability would require a set of engine-out 
sensors and the incorporation of gain scheduling. 

4. If it is deemed necessary to use more than the 5 
degrees of rudder needed by the RQAS, the actuator 
would need to be proportionately stronger. 

6.2 RECOMMENDED RESEARCH 

Further development of this project must include 
the following steps: 

1. Perform a preliminary failure analysis for 
rudder deflections greater than 5 degrees. 

2. Use the gains required at the most 
critical condition to see their effect in other 
conditions , 

3. Examine and evaluate options available to 
sense an engine-out, 

4. Investigate the possibility of increasing 
the amount of rudder deflection available to 
the RQAS, 

5. Design and build the necessary hardware, 

6. Perform a flight test of the system on the 
Cessna 402B. 
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APPENDIX A Cessna 402B Stability Derivatives 

This appendix gives the method used for estimating 
the airplane angle of attack and determining its 
stability and control derivatives. 

Because weight, ving area, and dynamic pressure are 
known, the airplane lift coefficient can be found from 
the equation: 

W - C L q S 

Prom this lift coefficient, the ACl due to flaps from 
reference 2 was subtracted. This lift coefficient was 
then found on Figure A.l and its corresponding angle of 
attack was read. The airplane angle of attack was then 
used to obtain the non-dimensional derivatives from 
reference 2. These are listed for an angle of attack, 
a, of 8.5° in Table A.l. 

The engine model was taken from reference 2. This 
gave the maximum power of a C402B engine as 300 brake 
horsepower. The propeller efficiencies, average and 
actual, are given in Figure A. 2. Thrust was then 
calculated using the following equation: 

t 

T a 550 BHP n p (A.l) 

V 

where 

np ■ Propeller Efficiency 
V ■ Airplane Speed 


BHP ■ Engine Brake Horsepower 





Table A.l Lateral-Directional Non-dimensional 
Stability Derivatives 


One Engine Out 
Pull Throttle 
Pull Flaps 


cy& 

a 

-0.670 

cm 

a 

0.195 

cy«DP 

a 

0.0 

Cn0 

a 

0.129 

Cn$R 

a 

-0.0795 

CnfiDP 

a 

0.0057 

ds 

a 

-0.0888 

C1 6R 

- 

0.0146 

C1 6DP 

- 

-0.0685 


Cessna 402B 
V = 130 fps 
Xjg a 0*25 

Derivatives 

rad” 1 

rad” 1 

rad” 1 

rad” 1 

rad" 1 

rad" 1 

rad" 1 

rad" 1 

rad" 1 


h * 0 ft 
W » 6200 lbs 


Cy p « -0.00063 
Cy r = 0.42 

Cn p • -0.084 
Cn r - -0.170 
Cl p « -0.81 
Cl r » 0.216 


V 


♦ 


A. 3 



Propeller Efficiency at Several Plight 
Conditions 



isumed Average Propeller Efficiency for All 
oi : rnnH i f inns 


Figure A.2 Propeller Efficiencies 
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From Figure 1, it can be seen that the thrust from one 
engine creates a moment about the airplane center of 
gravity with a moment arm of 7.50 feet. The orientation 
of the thrust line of one of the engines is given in 
reference 2. By finding the X and Z components of the 
thrust, the moments, H<p and L?, can be found from: 

N^> ® 7.50 T x , and 
L*r * 7.50 T z . 

Reference 2 also gives the AC^ due to an engine-out. 

ANd in Equation 2.1 is the yawing moment due to drag on 
the inoperative engine and is: 

AN d = 7.50 AC D q S 

Equation 2.1 was then solved to obtain the sideslip, 

differential flap deflection, and rudder deflection 

required. To account for changing the vertical tail 

size, the non-dimensional derivatives were recalculated 

and are summarized in Table A. 2. 

The dimensional derivatives were then calculated as 

shown in Table A. 3. These were calculated using: 

Ixx * 11100 slug ft 2 
♦ 

I 2 z » 14900 slug ft 2 
I xz * “583 slug ft 2 

To use these in the state matrices, Table A. 4, they must 
be in the form shown in Table A. 5, where, 

Ai “ ^xz^xx and ■ Ijcz^zz • 
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Table A. 2 Variation of Derivatives witn 
vertical Tail Size 


Cessna 402B (PCI) 


Relative 
V.T. Size 

S 

deg” 1 

C Vd F _ 

dec"* 

j S* 

i dec - 

! 

i dec"* 

0.50 

0.00085 

G.CC16S 

j - 0. 00070 

0 . 000155 

0.75 

0.00165 

0.00248 

i - 0.00104 

■ 0.000222 

1.00 

0.00240 

C.0023C 

| -0.00138 

0 . 000210 

1.25 

0.00216 

0.00413 

! -0.00174 

0.000288 


£? te : c no * -0.00062 deg~- was calculated f roir ref 
WB 
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Table A. 4 Lateral-Directional Dimensional Stability 

Derivatives [3] 


q l SC v 

1 y* -2 

Y a « ■■ (ft sec ) 

0 


q.SbC 

1 y „ -i 

Y_ » ~ (ft sec X ) 


p 2ml' 


q.SbC 

1 y r -1 

• “Suf <£t “ c > 

q.SC -2 

y o, (ft sec” 14 or 

Y m mmmuh* a»7 

6^ m ft sec deg ) 

h SC y -2 

y $, (ft sec 4 or 

v a ■ K _? -1 

6 m ft sec ~ deg ) 


q l SbC l 

L a * — (sec 4 ) 


“6 I 


xx 


q^b 2 ^ 

» -t; — rr- 2, (sec" 1 ) 


'P 21 U. 
v xx 1 


q 1 Sb 2 C i 

’ X77 (sec ' l) 


0l SbC tj 

5 a -*> 

L, » T (sec or 

* ** ..c- 2 4.,' 1 ) 


V bC t 


R , -2 

— (sec or 


\ \x -•> 

R xx sec * deg *) 


q SbC 

n a -? 

N a - * (sec ~) 


’a i 


22 


N - — (sec" 2 ) 

T 8 22 


<l Sb \ . 

N - (sec ) 


P 21 U. 
* 22 1 


q,Sb 2 C 

N - , 7— (sec l ) 


* 2I 22°1 


^l SbC n 




N 


a — - (sec or 

A 22 sec" 2 deg" 1 ) 


q,SbC 
1 n i 

6 R -2 

N. • — : (sec or 

R 22 sec" 2 deg" 1 ) 
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Table A. 5 Modified Lateral-Directional Dimensional 
Stability Derivatives [6] 


% 


b 


Yg » Yg/Ui 
Yp « Yp/Ui 
Yp » (Y r /Ui)-1 
Y$ = gcos8i/Ui 

Y *df * Y «df /Ul 

Y *sr a Y *sr /U * 

Lg » (AiNg+LgJ/d-AiBx) 
L p ' » (AiNp+LpJ/d-AiBi) 
I*r • (AlN r +L r ) /d-AiBj,) 


Lfi df - (AiM^^L^^/d-AxBx) 

W sr - (A 1 K« ir vLs ar )/(l-A 1 B 1 ) 

Ng » (BiLg+NgJ/d-AiBj.) 

V a (BiLp+NpJ/d-AiBi) 

Nf « (BiLr+NrJ/d-AiBx) 

N« df * 

N« S) . » (BlLg sr +Ng sr )/(l-A 1 B 1 ) 
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